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We investigate the T-S phase diagram of the electron-doped systems in high-T c cuprates where 
S is doping rate. Taking the ingap state and the super- exchange interaction J s , we calculate 
the superconducting transition temperature T c , the antiferromagnetic transition temperature 
Tn, the NMR relaxation rate 1/Ti with the antiferromagnetic fluctuations in the fluctuation- 
exchange (FLEX) approximation and with the superconducting fluctuations in the self-consistent 
i-matrix approximation. Obtained phase diagram has common features as those in the hole- 
doped systems, including the antiferromagnetic state, the superconducting state and the spin 
gap phenomenon. Doping-dependences of Tn, Tc and Tsg (spin gap temperature) are, however, 
different with those in the hole-doped systems. These differences are due to the intrinsic nature 
of the ingap states which are intimately related with the Zhang-Rice singlets in the hole-doped 
systems and are correlated d-electrons in the electron-doped systems, respectively, which has 
been shown in the d-p model. 

KEYWORDS: electron-doped cuprate, pseudogap, spin gap, superconductivity, superconducting fluctuation, spin 
fluctuation, superexchange interaction, d-p model 
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In the high-T c cuprate superconductors, the super- 
conducting (SC) states appear as holes or elections are 
doped into the antiferromagnetic (AF) insulatoiBcP . The 
SC states in the electron-doped cuprates have d p 
symmetry as those in the hole-doped cuprates dc 
which suggests existence of a common mechanism in 
both kinds of cuprates. The doping ranges in which 
the AF states and the SC states exist in the electron- 
doped cuprates are, however, defferent from those in 
the hole-doped cuprates. In the hole-doped cuprates, 
the AF states exist in narrow doping range (0 < S £ 
0.02), and the SC states exist in wide doping range 
(0.05 < S < 0.25). In the electron-doped cuprates, 
on the other hand, the AF states exist in wider range 
(—0.13 ;$ 5 < 0) and the SC states appear in narrower 
range (—0.18 % 5 ;$ —0.13) than those in the hole-doped 
cuprates. 

In our previous study we have obtained the T— S phase 
diagram in the hole-doped region with the AF fluctua- 
tions in the fluctuation-exchange (FLEX) approximation 
and the SC fluctuations in the self-consistent t-matrix 
approximations' , which have been based on the ingap 
state and the super-exchange interaction J s obtained in 
the d-p model by using tb.e_slave boson technique and 
the 1/A-expansion theoryu'ErEP . We have obtained the 
phase diagram including the AF state, the SC state with 
d a ,2_ v 2-symmetry and_th£_spin gap phenomenon induced 



by the SC fluctuation^ ' @ O , which are consistent with 
those observed in the hole-doped cupratesBlij). 

In the present study we investigate the T — 5 phase dia- 
gram in the electron-doped region by calculating the su- 
perconducting transition temperature T c , the antiferro- 
magnetic transition temperature Tn and \fT\T of NMR 
by using the same approximation as that taken in our 
previous studytf . Obtained phase diagram has common 
features as those in the hole-doped region, including the 
AF state, the SC state and the spin gap phenomenon. 
Doping-dependences of Tn, T c and Tsg (spin gap tem- 



perature) are, however, different from those in the hole- 
doped systems. The AF state exists approximately in 3 
times wider range than those obtained in the hole-doped 
region. The SC state, on the other hand, appears in 
narrower range. In addition we show that the spin gap 
phenomenon appears in narrower range than those ob- 
tained in the hole-doped region. The electron-hole asym- 
metry is due to the intrinsic nature of the in-gap state 
described in the followings, while the electron-hole asym- 
metry has been studied by single band models such as 
the t—J model or the single band Huh 
troducing the long range hoppingsOc 

We take the d-p model which is a kind of the two-band 
Hubbard model for describing the electronic system in 
the Cu02 plane: 
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which is treated within the physical subspace where local 
constraints Qi = J2a d ty^ + bfbi = 1 strictlx hold in 
order to exclude double occupancy of d-holesEr EHP . In 
the above, cj^, d^ and bi are annihilation operators for a 
p-hole with a wave vector k and spin <r, a pseudo fermion 
representing a single occupation of ri-hole of i-th site and 
a slave boson representing a vacancy of the ri-hole of i- 
th site, respectively, and t.jj = tj i exp(— ik • Rj), with 
ijj. = 2<d p [l — T;(cosk x a+cosk y a)]2 , where a is the lattice 
constant (we set a — 1), and Al is the total number of 
lattice sites. Here e v and Ed stand for the energies of the 
p- and ri-levels measured from the chemical potential /i, 
respectively, and td p stands for the d-p transfer. 

The 1/A-expansion is a suitable method for treat- 
ing correlations systematically. First, the strong lo- 
cal correlations are suitably included in the self-energy 
in the leading order of 1/A-expansion. The quasi- 
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particle Green's functions 
1/A-expansion are given by 



he leading order in the 



G {k,uj) = J2 A 7 (k)/(w-£ 7 (k) + iO H 

7=± 



(2) 



with £y(k) 



Abt 2 )V2] 



|[£p + UJQ + j((Sp - W ) 2 -T ^o k ; 

A 7 (k) = 7(£ 7 (k) ~ w )/(^+(k) - £L(k)), where 7 = - 
and + denote the ingap state and the p-band, and iV rep- 
resents the degeneracy of d-hole (JV = 2 in the present 
case). The binding energy ujq and the number b of 
the slave boson are self-consistcntly determined together 
with fi, and then /i is located in the ingap state. 

The system is the charge transfer (CT) type of the 
Mott insulator at_hole doping rate 5 = and in the 
region of A > AjE$ , where A = e p — Ed stands for the 
CT gap, and the total hole number is given by n = 1 + 5. 
In the hole-doped region, doped p-holes form the quasi- 
particle band called as the ingap state inside the CT 
gap near the p-band. It is due to the resonating state 
of the slave boson, which corresponds to the Zhang- Rice 
singletEfP. Those resuLts rJ ar£_ 1 consistent with the results 
of optical experimentsE^cJ'tZP . 

In the electron-doped region where vacancies of d-holes 
are doped, on the other hand, the ingap state is the d-like 
band inside the CT gap near the localized d-state. The 
band width and the intensity of the in-gap state approxi- 
mately equal loq and b, respectively, and are proportional 
to 1 6 1 when \S\ <C 1. The intensity in the electron-doped 
region is almost 3 times larger than that in the hole- 
doped region, i.e. the intrinsic electron-hole asymmetry 
induced by the local correlation on the d-sites. Those ba- 
sic features are not changed-by taking the higher order 
terms in the l/iV-expansionc3 ) . Those results are consis- 
tent with the results obtained recently by the dynami- 
cal mean field theory (DMFT) in the two-band Hubbard 
modelBBB. 

Next, concerning the intersite correlations, it has been 
shown that the super-exchange interaction J s is the most 
dominanLterm within 0(1/N) in the vicinity of the CT 
insulators' . We derive the coupled equations of the AF 
fluctuations and the SC fluctuations which are induced 
in the ingap state by J s . The AF fluctuations via J s are 
given by 



V(q,w) = J s (q)/(1 ~ J s (q)x S (q,^)), 



(3) 



with J s (q) = — J s (cos(q x ) + cos(q y )) and the irreducible 
spin susceptibility 



X 



'(q,w) =^L 1 E fi k + qk^ 1 / dx f( x ) 
k 

x [ImG(k + q, x)G{k, -w + x)* 
+ G(k + q, w + x)ImG(k, x)], 



(4) 



where tt k k , = fe 2 ^,/((i?_ (k) - u f{E^k') - w ) 2 ) 
and f(E) = (exp{£/T} + l)" 1 . 

It was shown that a component with the d x 2_ y 2 sym- 
metry among various components of the spin-fluctuation- 
mediated interaction contributes dominantly to the pair- 
ing interaction^ We take only the component with the 



d x 2_ y 2 symmetry of the pairing interaction as 

Vd 



A" 2 



(5) 
k k' 

with ipjf = cos(k x a) — cos(k y a). 

The irreducible pairing susceptibility x p (q, w) is given 
by 



X 



>(q, W ) = -At 1 J2 ViAk ^kvr- 1 J dx 



x [/(i)ImG(k, x)G(q -k,ui- x) 
- f(-x)G(k, oj - x)ImG(q - k, x)] 



(6) 



and A k k , = b 2 t?t?/((irT) 2 + luq) 2 . In the above equa- 
tions, both ^ujj/ and A k k , are the vertices connecting 
p- and d-holesB 1 

Self-energy corrections due to the SC fluctuations in 
the i-matrix approximation are given by 

Escf (k, w) = iffcJVf 1 J2 A kq k^ 1 / dx 



x [/(ar)T(q, x + w)ImG(q - k, x) 
- fl (x)ImT(q,x)G(q-k,x-w)% 
with the i-matrix of the SC fluctuations 

T( q , lj) = ^x p (q, w)/(i - «dx p (q, w)), 



(7) 



(8) 



where g(E) = (exp{E/T} - l)" 1 . 

Self-energy corrections due to the AF fluctuations 
in the fluctuation exchange (FLEX) approximation are 
given by 

S APf (k, w) =N£ 1 Y J A k.k q^ * / dx 

q J 

x [g(x)lmV(q,, x)G(k — q, uj — x) 

-f(-x)V(q,uj- x)ImG(k - q,x)]. (9) 

The renormalized Green's function is given by 

G(k,w) = [Go(k^)- 1 

-S SCf (k,w)-E AF£ (k,a;)]- 1 , (10) 

Both the SC and AF fluctuations are self-consistently 
treated by solving the coupled equations (3)-(10). 

Using the solutions of the above coupled equations, we 
calculate the NMR relaxation rate, 

l/T.T = AT- 1 VF ak (q) lim , Im X s (q, u)/u, (11) 

* — * u— >0 

g 

where x s (q,^) = J s (q)x s (q, w)) andF a6 (q) 

are the form factorsH. In the cuprates, \jT\T indicates 
the degree of the AF fluctuation. 

We define T c as the temperature determined by the 
condition, 



1-^(0,0) -a = 



(12) 



, where a represents the degree of the 3-dimension (3D) 
effect in the quasi-2D electronic system. In the present 
model, 1 — VdX p {0, 0) does not reach zero at finite tern- 
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peratures by the effect of Escf (k, w) , because we take a 
2D electronic system in the calculation of the above cou- 
pled equations. We take a = 0.5 in the present study so 
as to have finite T c of O(100K). 

We define Tn as the temperature determined by the 
condition, 

C = l- J s (Q)x s (Q,0)-/3 = (13) 

, where Q = (tt, tt) and (3 represents the degree of the 
3D effect. In the present model 1 — J s (Q)x s (Q, 0) also 
does not reach zero at finite temperatures by the effect of 
SAFf (k, oj) in the 2D electronic ststem. We take (3 — 0.02 
in the present study. 

We define To as the temperature at which the AF 
fluctuations begin to be appreciable (^(1/TiT)|t=t = 
—50(1/TiT)|t=o. 08, 5=o.i in the present study). It is as- 
sumed that To is the upper boundary ofihe .anomalous 
metallic phase in the hole-doped cuprateoli3'E-$ . 

In actual numerical calculations, throughout the 
present study, we set 2td p = 1.0 (eV), A = 2.5 and 
J s = 0.1. For example, we have ojo — 0.0780 and 
b = 0.0515 at 6 — 0.1. The total number of discrete 
points taken for the q-summation over the 2D first Bril- 
louin zone is 32 x 32. The w-integral over the region 
from — 2dj to 2o>o is replaced by the w-summation of 80 
discrete points. 

Figure 1 shows the T-5 phase diagram obtained by 
solving the above coupled equations, where the phase 
diagram in the hole-doped region (5 > 0) has been 
obtained in our previous study©. The AF state in 
the electron-doped region persists to higher doping rate 
S = —0.13 than the doping rate S = 0.05 where that 
in the hole-doped region persists. The SC state appears 
in narrower doping range —0.19 ,$ 8 ^ —0.13 than the 
doping range where that in the hole-doped region ap- 
pears. In addition it is shown that the spin gap ap- 
pears in the doping range —0.16 ^ 5 ^ —0.13, which 
is consistent with the results observed by NMR very re- 
ccntlyEP . Those features account for the phase diagrams 
observed both in the electron-doped cuprates and in the 
hole-doped cuprates. 

The AF fluctuations dominate in the region of Tsg ~ 
T < T , while the SC fluctuations dominate in the re- 
gion of T c < T < TsqEP . This crossover behavior induces 
the anomalous T-dependence of 1/T\T as shown in Fig. 
2. As a result of the domination by the SC fluctuations, 
T c has a maximum at 5 = 0.1 and decreases with S de- 
creasing in the hole-doped region. If the AF state is re- 
moved, T c in the electron-doped region has a maximum 
at S = —0.1 and decreases as \5\ decreases. 
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Fig. 1. The T-8 phase diagram where the system is in the super- 
conducting state with d x 2_ y 2 symmetry in the region of T < T c , 
and is in the antiferromagnetic state in the region of T < Tjq. 
The spin gap temperature Tsg is denned as the temperature at 
which 1/TiT has a maximum as shown in Fig. 2. The metallic 
region in T To is the anomalous metallic phase. 
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Fig. 2. T-dependences of 1/TiT at <5 = -0.08, -0.11, -0.12, 
-0.13, -0.14, -0.15, -0.16, -0.18 and -0.2. The triangles 
represent the spin gap temperature TgQ at which l/T\T has a 
maximum. The solid arrows and the dotted arrows represent T c 
and Tjm, respectively. 



Figure 2 shows T-dependences of 1/TiT. At 6 > 
—0.12, 1/TiT increase as T decreases until the system 

reaches T = T N . At S = -0.13 0.16, l/T x T has 

a maximum at the spin gap temperature TgQ and de- 
creases as T decreases in the region of T c < T < Tsg, 
which is the spin gap phenomenon induced by the SC 
fluctuation. If the AF state were to be removed, \jT{T 
has a maximum at Tsg m the region of 5 > —0.12, and 
then Tsg continues to increase with \5\ decreases. 



Figure 3 shows T-dependences of r defined in eq. (12) 
(the solid lines) and £ defined in eq. (13) (the dotted 
lines). At 5 — —0.13 ~ —0.18, r decreases as T decreases 
and reaches zero at T — T c . At 5 £ — 0.12, £ decreases as 
T decreases and reaches zero at T = Tn, while £ does not 
reach zero in the region of 6 < —0.13. The upturn at S = 
—0.13 and the saturation at 6 = —0.14 ~ —0.16 in the T- 
dependences of £ are mainly due to the SC fluctuations, 
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while those at 6 ^ —0.16 are due to the decrease of the 
AF fluctuation. 




Fig. 3. The solid lines represent T-dependences of r at <5 = —0.2, 
-0.18, -0.16, -0.15, -0.14 and -0.13, which reach zero at T = 
T c . The dotted lines represent T-dependences of £ at <5 = —0.2, 
-0.18, -0.16, -0.15, -0.14, -0.13, -0.12, -0.11 and -0.08, 
which reach zero at T = Tjm. 



In the phase diagram obtained in the present study, T c 
and Tn depend on the phcnomenological parameters a 
and (3, respectively, representing the degree of the 3D ef- 
fect depending materials. The qualitative features of the 
(5-dependences of T c and T N are, however, independent 
of the choice of these parameters, and TgQ is determined 
without indefiniteness. By solving the coupled equations 
(3)-(10) for the anisotropic 3D system, ambiguities due 
to these parameters should be excluded. 

In summary, we have obtained a unified phase diagram 
accounting for the phase diagrams observed both in the 
electron-doped cuprates and in the hole-doped cuprates 
including the AF states, the SC states and the spin gap 
phenomena. The AF fluctuations and the SC fluctua- 
tions induced in the ingap states by J s are the key fac- 
tors for describing those, especially the spin gap phenom- 
ena. The electron-hole asymmetry on the 5-dependences 
of Tn, T c and Tgc, is the result from the differences of 
intrinsic nature of the ingap states due to the strong cor- 
relation on the d-sites in the d-p model. 

The present work has been partially supported by the 
Grant-in-Aid for Scientific Research from the Ministry 
of Education, Culture, Sports, Science and Technology, 
Japan. 



[6] A. Kobayashi, A. Tsuruta. T. Matsuura and Y. Kuroda: J. 

Phys. Soc. Jpn. 70 (2001) 1214. 
[7] H. Jichu, T. Matsuura and Y. Kuroda: J. Phys. Soc. Jpn. 

59 (1990) 2820. 

[8] Y. Ono, T. Matsuura and Y. Kuroda: J. Phys. Soc. Jpn. 64 
(1995) 1595. 

[9] S. Fukagawa, A. Kobayashi, K. Miura, T. Matsuura and 

Y.Kuroda: J. Phys. Soc. Jpn. 67 (1998) 3536. 
[10] A. Kobayashi, A. Tsuruta. T. Matsuura and Y. Kuroda: J. 

Phys. Soc. Jpn. 68 (1999) 2506. 
[11] A. Kobayashi, A. Tsuruta. T. Matsuura and Y. Kuroda: J. 

Phys. Chcm. Sol. 62 (2001) 199. 
[12] Y. Kosuge, H. Sato, A. Kobayashi, T. Matsuura and Y. 

Kuroda: J. Phys. Chcm. Sol. 62 (2001) 243. 
[13] G.-q. Zheng, H. Ozaki, W. G. Clark, Y. Kitaoka, P. Kuhns, 

A. P. Reyes, W. G. Moulton, T. Kondo, Y. Shimakawa and 

Y. Kubo: Phys. Rev. Lett. 85 (2000) 405. 
[14] Y. Suzumura, Y. Hascgawa and H. Fukuyama: J. Phys. Soc. 

Jpn. 57 (1988) 2768. 
[15] T. Tohyama and S. Mackawa: Phys. Rev. B 49 (1994) 3596. 
[16] K. Kuroki and H. Aoki: J. Phys. Soc. Jpn. 67 (1998) 1533. 
[17] H. Kontani, K. Kanki and K. Ucda: Phys. Rev. B 59 (1999) 

14723. 

[18] H. Kondo and T. Moriya: J. Phys. Soc. Jpn. 68 (1999) 3170. 
[19] D. Manske, I. Eremin and K. H. Bennemann: Phys. Rev. B 

62 (2000) 13922. 
[20] Y. Yanasc and K. Yamada: J. Phys. Soc. Jpn. 70 (2001) 

1659. 

[21] P. Coleman: Phys. Rev. B29 (1984) 3035. 

[22] B. Jin and Y. Kuroda:_ J. Phys. Soc. Jpn. 57 (1988) 1687. 

[23] D. S. Hirashima, Y. Ono, K. Miura, T. Matsuura and Y. 

Kuroda: J. Phys. Soc. Jpn. 60 (1991) 2269. 
[24] F. C. Zhang and T. M. Rice: Phys. Rev. B 37 (1988) 3759. 
[25] S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura and S. 

Tajima: Phys. Rev. B 43 (1991) 7942. 
[26] C. T. Chen, F. Sette, Y. Ma, M. S. Hybertsen, E. B. Stechel, 

W. M. C. Foulkcs, M. Schultcr, S. W. Chcong, A. S. Cooper, 

L. W. Rupp, Jr., B. Batlogg, Y. L. Soo, Z. H. Ming, A. Krol 

and Y. H. Kao: Phys. Rev. Lett. 66 (1991) 104. 
[27] N. B. Brookes, G. Ghiringhelli, O. Tjernberg, L. H. Tjeng, 

T. Mizokawa, T. W. Li and A. A. Mcnovsky: Phys. Rev. Lett. 

87 (2001) 237003. 
[28] A. Tsuruta, A. Kobayashi, T. Matsuura and Y. Kuroda: J. 

Phys. Soc. Jpn. 69 (2000) 3342. 
[29] T. Mutou and D. S. Hirashima: J. Phys. Soc. Jpn. 66 (1997) 

2781. 

[30] Y. Ohashi and Y. Ono: J. Phys. Soc. Jpn. 70 (2001) 2989. 

[31] Y. Ohashi and Y. Ono: private communication. 

[32] M. Azami, T. Matsuura and Y. Kuroda: J. Phys. Soc. Jpn. 

66 (1997) 2811. 
[33] T. Imai: J. Phys. Soc. Jpn. 59 (1990) 2508. 
[34] T. Nishikawa, J. Takeda and M. Sato: J. Phys. Soc. Jpn. 63 

(1994) 1441. 

[35] Y. Kitaoka and G. -q. Zheng: private communication. 



[1] M. Sato: Physica C 263 (1996) 271. 

[2] H. Takagi, Y. Tokura and S. Uchida: Physica C 162-164 
(1989) 1001. 

[3] C. C. Tsuei and J. R. Kirtley: Phys. Rev. Lett. 85 (2000) 
182. 

[4] N. P. Armitage, D. H. Lu, D. L. Feng, C. Kim, A. Damascelli, 
K. M. Shen, F. Ronning, Z. -X. Shen, Y. Onose, Y. Taguchi 
and Y. Tokura: Phys. Rev. Lett. (2001) 86 1126. 

[5] T. Sato, T. Kamiyama, T. Takahashi, K. Kurahashi and K. 
Yamada: Science 291 (2001) 1517. 



